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Abstract 

In gapped spin ladder and spin-Peierls systems the introduction of 
disorder, for example by doping, leads to the appearance of low energy 
midgap states. The fact that these strongly correlated systems can be 
mapped onto one dimensional noninteracting fermions provides a rare 
opportunity to explore systems which have both strong interactions 
and disorder. In this paper we will show that the statistics of the zero 
energy midgap wave functions V'o(^) in these models can be effectively 
described by Liouville Quantum Mechanics. This enables us to calcu- 
late averages over disorder of the products iI)q{xi)iI)q{x2)...'4'q{xn) (the 
explicit calculation is performed for N = 2,3). We find that whilst 
these midgap states are typically weakly correlated, their disorder av- 
eraged correlations are power law. This discrepancy arises because the 
correlations are not self-averaging and averages of the wave functions 
are dominated by anomalously strongly correlated configurations; a 
fact which is not always appreciated in the literature. 
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1 Introduction 



One dimensional quantum spin systems have long fascinated both theorists 
and experimentalists; more recently there has been particular interest in 
the behaviour of such systems in the presence of disorder (see eg. and 
references therein). There are several reasons for this, but one significant mo- 
tivation is the comparatively recent experimental realisation of doped quasi- 
one-dimensional spin-Peierls and spin ladder systems 0. One of the most 
important general results in these materials is the appearance upon doping 
of magnetic states at energies well below the spin gap- which dramatically 
changes their magnetic properties (see eg. |^ and references therein). 

As well as the new results emerging from experiment, theoretical work 
has revealed many unusual features in a range of different one dimensional 
disordered spin models. Prominent among these results is the occurrence of 
Griffiths phases; here the low energy response is dominated by anomalously 
strongly correlated regions of the system. The property related to that is a 
sharp distinction between typical and disorder averaged correlations; the lat- 
ter are much stronger and typically power law at criticality . The fact that 
the correlation functions are not "self-averaging" means that rare strongly 
correlated configurations dominate the disorder averaged quantities. Thus 
the disorder averaged correlations can be much stronger than one would 
naively guess based on the implicit assumption that "averaged" and "typ- 
ical" are the same thing. This is important because it is these averaged 
correlations which will be relevant in experiments. 

Our approach to these one dimensional spin systems is to exploit the 
various mappings to fermionic theories. This provides an alternative way of 
thinking about the problem which can be related to the extensive body of lit- 
erature on disordered fermion systems. The spin sector of both the two chain 
spin ladder ||^ and the spin-Peierls system |^ can be represented in terms of 
massive noninteracting fermions. When these systems are lightly doped their 
conductivity shows variable range hopping behaviour, which means that the 
charge carriers are strongly localised and no band forms. Thus to investigate 
magnetic properties it is reasonable to continue to work with a pure spin 
model without fiuctuations in the charge sector. 

We will argue that doping with nonmagnetic impurities can be repre- 
sented by a type of randomness in the mass of the fermions of the effective 
theories 0,0. Specifically the mass has fixed magnitude mo (the value of 
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the undoped spin gap) except that it flips between +mo and —mo at the lo- 
cations of impurities such that (m) = 0. This nontrivial background admits 
localised low energy states well below the gap in addition to the massive bulk 
modes. 

It is useful to comment briefly on some insights from the theory of elec- 
trons in disordered metals, in order to better understand the connection with 
spin theories. One of the reasons that this problem is so complex is the ex- 
istence of many different length scales. A great deal is known about the 
behaviour of metallic-type extended states which explore the whole sample, 
i.e. the limit of small wavefunction amplitudes t = \iIj{x)\'^; this depends only 
upon the global symmetry of the ensemble and can be derived from random 
matrix theory (see eg. ^ and references therein). 

However, at length scales much less than the localisation length but 
much larger than the mean free path /, strongly localised states with anoma- 
lously high local amplitude are important (in the conventional localisation 
this situation emerges when disorder comes from scattering on magnetic im- 
purities). These states look far from metallic in the region I <^ L < Lc and 
their correlations are sensitive to local variations in the potential; this is es- 
sentially the origin of the distinction between typical and disorder averaged 
correlation functions in this regime 0. This is in sharp contrast to the be- 
haviour of metallic states which are equally affected by the random potential 
at all points in the sample, and not surprisingly the wavefunction statistics 
are no longer simple. Here we can see a relation with the behaviour of Id 
disordered spin systems; it is this kind of state which corresponds to localised 
low energy spin degrees of freedom. 

Recent work on this problem using a variety of approaches has shown 
that the wavefunction statistics in this regime can be described by Liouville 
field theory |p and that, for example in the two dimensional case relevant to 
the quantum Hall effect, one observes such interesting phenomena as multi- 
fractality 

It was pointed out in that in a system of one dimensional fermions with 
random mass m and (m) = 0, the Lyapunov exponent 7(£^) ~ l/(— Inii^) 
tends to zero as E ^ and so for low energies we are always at length 
scales much less than the localisation length Lc ~ I/7. Thus it is not at all 
surprising that in one dimensional systems of length L > I we see a departure 
from the universal features predicted by random matrix theory 0. 

In this paper we will study one dimensional Dirac fermions with a random 
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mass. We show that the correlation functions of the so-called prelocalised 
zero energy states can be calculated using as an effective theory Liouville 
quantum mechanics. This analysis is rather general and does not depend 
upon the specific form of the disorder. In section (3) we will then describe 
the application of these results to spin systems; here the low energy states 
represent the midgap magnetic states mentioned above. The advantage of 
this approach to the midgap states is that it enables us to get a clear picture 
of their behaviour and specifically their correlations, information which is 
hard to obtain using more traditional methods. 



2 Effective Theory 

Suppose we have a system of 1+1 dimensional non- interacting Dirac fermions 
with a random position dependent mass. The Hamiltonian can be written; 

H = I dxip^ {-ta^-^ + a^m{x)j^ (1) 





where and are the standard definitions of the Pauli matrices. This 
leads to the equations of motion, in component form; 

— I" "^(^)^ '"(^) = Ev{x) 

^ v{x) = Eu{x) (3) 

Even if m{x) is finite and nonzero almost everywhere, if it passes through 
zero at some points there will be normaliseable bound states with energy 
close to zero- midgap states. This can be seen from replica, supersymmetry 
and numerical calculations where a peak appears in the density of states at 
E = (see eg. |^ and references therein). However, within these approaches 
it is notoriously difficult to calculate the correlation functions. 

It is remarkable that one of the low energy solutions of eqs.(^ can be 
explicitly written down; it is the zero energy eigenmode; 

Mx) = J )exp(-\/(x)) (4) 
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V{x) = / m{y)dy (5) 







L 



= / dxexp{-2V{x)) (6) 







Here it is our purpose to calculate correlations of these zero energy eigenstates 
for the case where (m), the average mass is zero. 

In disordered systems, we are usually interested in disorder averaged cor- 
relations. Here as always the difficulty is to correctly account for the normal- 
isation factor (denominator) in the average. Following Kogan, Mudry and 
Tsvelik p we introduce a dummy variable /i which allows us to transfer all 
the V{x) dependence from the denominator into the numerator; a disorder 
average can then be performed. 

In the continuum limit, it is apparent that the quantity V{x) will 
behave like the position variable in a random walk where x is the "time" 
coordinate. Therefore, going to a path integral formulation we find in the 
continuum limit the following Gaussian probability distribution for V (for 
more details see lITOi); 



P{{V})DV = exp 



dx{d,V) 



DV (7) 



1 

2^ 

where g parametrises the strength of the disorder. This will be valid for 
sufficiently long distances in the correlation functions; different realisations 
of disorder may give rise to different behaviours at small separations. In 
the case of "telegraph disorder" where m randomly assumes either ±mo, 
9 = — where oq is the lattice spacing and v is the velocity in the Dirac 
equation (^ (often taken as 1 for simplicity). For more general forms of 
disorder, g will roughly speaking be larger for broader distributions of m. 
Note that g is independent of impurity density- this will however affect the 
length scales above which the distribution (^ is valid. 

Then the disorder averaged correlator {{X) denotes the disorder average 
of the quantity X) ; 



{\M^,)\\..\M^N)n = I DVP{{v})\M^,)\' \Mxn)\' (8) 

can be expressed as 



1=1 
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(9) 



where a is introduced to regularise the [i integral and 



dx 



•^9 



(10) 



is the action of Liouville quantum mechanics [|13[ where x is the "time" 
coordinate. 

We can now go to the canonical Hamiltonian form; 



H 
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2dV'' 



+ /xexp(-2V) 



(11) 



and find the eigenstates of the effective theory, which allows us to calculate 
the correlation functions. The detailed calculation of the one, two and three 
point correlators is described in the Appendix. In fact all of the correlators 
are proportional to a constant factor which depends upon an integral over 
the dummy variable /i (which must be suitably regularised) and an unspec- 
ified normalisation factor A. Since the same factor appears in front of all 
the correlators, we may fix its value by insisting that the wavefunctions are 
correctly normalised: 



hence 



dx{\ij,{x)\^) 



(l^o(O)H 



1 
L 



(12) 



(13) 



where L is the size of the system. Comparing this with the result of A2, Eq. 
( |33D we find for the two and three point correlators the following correctly 
normalised power law behaviour: 



(|^o(x2)|Vo(^2)r) 
{\M^z)?\M^2)?\U^l)?) 




(14) 

(15) 
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where we have chosen without loss of generahty X3 > X2 > xi and xtj = 
Xi — Xj, and the derivation assumed that \gxij\ ^ 1, i 7^ j. This is a 
very clear example of the important difference between typical and disorder 
averaged correlations. From the expression (|D we would expect the typical 
correlations to decay exponentially. What happens is that atypical strongly 
correlated configurations dominate the average and give much stronger power 
law correlations in the disorder averaged quantities. 



3 Application to Spin Systems 

We discussed in the introduction some of the characteristic features of one 
dimensional disordered spin systems including the occurrence of Griffiths 
phases, and the related distinction between typical and averaged correlation 
functions. As we have seen a similar picture holds for the correlations of low 
energy fermion modes in one dimension; and this is not a coincidence but 
an expression of the deep relationship between certain one dimensional spin 
systems and systems of fermions. 

Since the analysis above is equally valid for Majorana (real) fermions, one 
can first of all make the direct observation that this model with random mass 
is equivalent to the quantum Ising model with random bonds. This has been 



studied by many authors (including eg [0]). In particular Shankar and 
Murthy fl^ found that the typical correlation functions of fermion bilinears 
decayed exponentially as the square root of distance as one would expect; our 
result shows in addition that at criticality the disorder averaged correlation 
functions of zero energy modes are power law. 

Let us now consider the doped spin-Peierls and spin ladder systems. Fab- 
rizio and Melin [^] have suggested that one can gain insight into the spin- 
Peierls system by considering the XY version- this maps to to a system of 
noninteracting massive Dirac fermions whose mass is proportional to the 
dimerisation. They argue that the introduction of impurities can be mod- 
elled simply by the presence of domain walls between vacua with dimerisation 
of opposite sign. This leads to a model of fermions with a mass which flips 
between +mo and —mo at impurities (by symmetry (m) = 0); one finds low 
energy midgap states localised in regions with impurities, of which the zero 
energy eigenstate considered in this paper is one example. 

In the language of spins these midgap states correspond to the effective 



7 



spins which appear at breaks or discontinuities in gapped spin chains [0], 
since for this model the z component of the magnetisation is proportional 
to the fermion density. 

These impurity spins are often treated in the literature as though they 
are essentially uncorrelated, based on the argument that they are typically 
weakly coupled for low impurity density. We emphasise again that whilst 
this is true, our work shows that their disorder averaged correlations are 
power law; and these are the relevant quantities in experiment. It is for this 
reason, for example that they do not give a simple Curie contribution to the 
susceptibility; Fabrizio and Melin showed that xiT) ~ l/(Tln^ T) a result 
which differs significantly at low temperatures from the noninteracting 1/T 
result. 

Of course in reality the problem of spin-Peierls systems is more complicated- 
there are very important three dimensional effects and thermal fluctuations of 
the lattice distortion. It is also unclear, even in an inherently one dimensional 
system, whether the introduction of impurities can really be represented so 
simply . Nonetheless we would argue that we have captured a crucial fea- 
ture; which is that at boundaries between the two degenerate ground states 
(domain walls = kinks in the fermion mass) we find localised magnetic low 
energy degrees of freedom. 

Similar arguments applied to a fermionic model of the two chain spin 
ladder where doping is introduced in the form of static kinks in the charge 
field, leads again to a model of fermions with mass fiipping between —mo 
and mo at the location of impurities. This can be seen from the bosonisation 
formalism, treating nonmagnetic impurities in a way suggested by Fukuyama 
et al. [^. At half filling in a spin ladder system, the charge modes will be 
frozen (gapped). Nonetheless, because the charge density is proportional to 
the gradient of the charge field, a nonmagnetic charged impurity (hole) can 
be represented as a static kink in the field. This assumes that there are 
sufficiently few impurities and that they are sufficiently strongly localised 
so that no band forms- these conditions appear to be satisfied in typical 
experimental systems 0. 

By taking these kinks into account in the effective theory of [^], it can 
be seen that their effect is to introduce kinks in the mass of the fermions 



describing the spin sector 12 



The interpretation in terms of spins is much as for the spin-Peierls sys- 
tem. In the effective theory M, the slow component of the magnetisation is 
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essentially proportional to the fermion density, and so these localised states 
can again be interpreted as low energy spin degrees of freedom. 

4 Conclusions 

We have shown that Liouville quantum mechanics can usefully be applied 
to calculate the properties of the zero energy localised states in a model 
of Dirac fermions with random mass. We have emphasised the relation of 
this model to some one dimensional disordered spin systems; one of the most 
important points is that at the critical point ((m) = 0) we find power law cor- 
relations for the disorder averaged correlation functions, even though typical 
correlations decay exponentially. This shows that even if midgap states are 
typically weakly correlated, anomalously strongly correlated configurations 
can dominate in physical quantities. 

Finally, it is very interesting in itself that one can derive an effective 
theory for midgap states with nontrivial, calculable correlation functions; 
although a limitation of our approach is that it only gives information about 
states with energy E = 0. It would be very interesting if the current approach 
could be extended to investigate states with finite energy. 
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A Calculation of Correlators 
A.l Two Point Correlator 



C21 



{\Mx2)\'\M^i)f) 



(16) 
(17) 
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Mo 



OO /■ + 00 



— oo J — OO 



dE 



+00 



dV{^\e 



U) 



-Ex2i 



(18) 
(19) 



where E = gv? /2, we have chosen without loss of generahty X21 = X2 — X1 > 
and the energy normahsed eigenstates of Liouville quantum mechanics at 
X = are given as follows [|^; 



|0) 



\u) 



AKo I W=^e-^ 



/2/i 
9 



27ig 



sinh TTU 




(20) 
(21) 



where A is a normalisation which we can fix later. The "time" (x) evolution 
of these states is denoted by subscripts on the bras and kets: 



-Exi 



u{V = V^) 



(22) 



First let us work out the integral over V, making the substitution x 
/W^)e-^; 



/ 
/ 



dV{0\e-^\u) 



00 
Ag2_ 



[sinh TTU 



xKo{x)Kiu{x)dx 



(23) 
(24) 



Which gives ||T5[| ; 



Ag 



[sinh 7TU 



TTU 
2 



.sinh^^ 



So now 



(25) 



Mo 



57r/i^ Jo 



udu sinh ttm 



TTU 

2 



sinh^^ 



exp( 



9X21 2^ 

— —u^i 



(26) 



For long distances gx2i ^ 1 the Gaussian factor in the integrand is very 
rapidly decaying, so the other parts of the integrand can be expanded for 
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small u, giving; 



j4V 
647r//2 



f 

Jo 



u exp( —u )du 



^ 3 
1 \2 



64V27r//2 \gx2i) 
which gives for the two point correlation function at long dist 



2„i 



21 



64V27r 



Jo u 



-afi 



li 



1 \i 



X2l 



A. 2 Normalisation 

To set the normalisation let us consider 



where 



Ml 



Ci = (iV'o(o)r) 

/ + 00 
dV^(0|e-2^|0) 
-oo 



+ 00 



-2y 



Again making the substitution x = {J2jji/g)e ^ we find that 



(iv^o(o)r 



f 
Jo 



— ( 
1^ 



-afi 



A. 3 Three Point Correlator 



C-. 



321 



{\M^3)\'\M^2)\'\M^i)\' 



C321 = / d/i Ms 
Jo 2 
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Inserting the resolution of unity as before; 



M3 
h 
h 
h 



dE^dE^ihhh 



00 

+cx> 



—00 
00 



dVi{Q\e 
dV2{ui\e 
dV'i{u2\e 



-2y2(0)|^_^^g(Sl-S2)x2 



|0)e 



-B2X3 



(36) 
(37) 
(38) 
(39) 



where Ei = {g/2)uf, i = 1,2. The integrals Ji and J3 are the same as 
eq. ([23|) whereas I2 is a little different- introducing again the variable x = 
/V^)e-^; 



1 \ 1 r+<^ 

— e(^i--^2)'^2(sinh7rMisinh7ru2)2 / xdxKiu^{x)Kiu^{x) (40) 

KU Jo 

2 2 

e^''^-''^^^^ (sinh nui sinh 7^2)'^ ~ ""^'^^ ^ (41) 



Ann 

iTT 



32/i ^ " cosh ui — cosh M2 

Choosing without loss of generality X21 = X2 — xi > and X32 = 2:3 — 0:2 > 
we can use the same argument as before, that for large separations the 
exponentials of energy will decay very rapidly, and so the other parts of the 
integrand (|36| ) can be expanded for small Ui,U2- Then we find; 



2„3 



M3 



102477/^3 Jo 



u^dui exp( 



+00 



u\du2 exp( 



gX32U2 



142) 



And we obtain for the 3 point correlation function at long distances; 



G 



A' 



321 



4096 



dfi 
— ( 



-ail 



\X21J \Xs2/ 



(43) 
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